Gammaretroviruses, novel viruses and pathogenic bacteria in Australian bats with neurological signs, pneumonia and skin lesions.

More than 70 bat species are found in mainland Australia. While most studies of bat viromes focus on sampling seemingly healthy individuals, little is known about the viruses and bacteria associated with diseased bats. We performed traditional diagnostic techniques and metatranscriptomic sequencing on tissue samples from 43 Australian bats, comprising three flying fox (Pteropodidae) and two microbat species experiencing a range of disease syndromes, including mass mortality, neurological signs, pneumonia and skin lesions. Of note, we identified the recently discovered Hervey pteropid gammaretrovirus in a bat with lymphoid leukemia, with evidence of replication consistent with an exogenous virus. The possible association of Hervey pteropid gammaretrovirus with lymphoid leukemia clearly merits additional investigation. One novel picornavirus and at least three new astroviruses and bat pegiviruses were also identified in a variety of tissue types, as well as a number of likely bacterial pathogens or opportunistic infections, most notably Pseudomonas aeruginosa.

Substantial viral and bacterial diversity at the bat-tick interface.

Ticks harbour a high diversity of viruses, bacteria and protozoa. The soft tick Carios vespertilionis (Argasidae) is a common ectoparasite of bats in the Palearctic region and is suspected to be vector and reservoir of viruses and other microbial species in bat populations, some of which may act as zoonotic agents for human disease. The Soprano pipistrelle (Pipistrellus pygmaeus, Vespertilionidae) is widely distributed in Europe, where it can be found inside or close to human habitation. We used meta-transcriptomic sequencing to determine the RNA virome and common microbiota in blood-fed C. vespertilionis ticks collected from a Soprano pipistrelle bat roosting site in south-central Sweden. Our analyses identified 16 viruses from 11 virus families, of which 15 viruses were novel. For the first time in Sweden we identified Issuk-Kul virus, a zoonotic arthropod-borne virus previously associated with outbreaks of acute febrile illness in humans. Probable bat-associated and tick-borne viruses were classified within the families Nairoviridae, Caliciviridae and Hepeviridae, while other invertebrate-associated viruses included members of the Dicistroviridae, Iflaviridae, Nodaviridae, Partitiviridae, Permutotetraviridae, Polycipiviridae and Solemoviridae. Similarly, we found abundant bacteria in C. vespertilionis, including genera with known tick-borne bacteria, such as Coxiella spp. and Rickettsia spp. These findings demonstrate the remarkable diversity of RNA viruses and bacteria present in C. vespertilionis and highlight the importance of bat-associated ectoparasite surveillance as an effective and non-invasive means to track viruses and bacteria circulating in bats and ticks.

Meta-transcriptomics reveals potential virus transfer between Aedes communis mosquitoes and their parasitic water mites.

Arthropods harbor a largely undocumented diversity of RNA viruses. Some arthropods, like mosquitoes, can transmit viruses to vertebrates but are themselves parasitized by other arthropod species, such as mites. Very little is known about the viruses of these ectoparasites and how they move through the host-parasite relationship. To address this, we determined the virome of both mosquitoes and the mites that feed on them. The mosquito Aedes communis is an abundant and widely distributed species in Sweden, in northern Europe. These dipterans are commonly parasitized by water mite larvae (Trombidiformes: Mideopsidae) that are hypothesized to impose negative selection pressures on the mosquito by reducing fitness. In turn, viruses are dual-host agents in the mosquito-mite interaction. We determined the RNA virus diversity of mite-free and mite-detached mosquitoes, as well as their parasitic mites, using meta-transcriptomic sequencing. Our results revealed an extensive RNA virus diversity in both mites and mosquitoes, including thirty-seven putative novel RNA viruses that cover a wide taxonomic range. Notably, a high proportion of viruses (20/37) were shared between mites and mosquitoes, while a limited number of viruses were present in a single host. Comparisons of virus composition and abundance suggest potential virus transfer between mosquitoes and mites during their symbiotic interaction. These findings shed light on virome diversity and ecology in the context of arthropod host-parasite-virus relationships.

Faecal virome of the Australian grey-headed flying fox from urban/suburban environments contains novel coronaviruses, retroviruses and sapoviruses.

Bats are important reservoirs for viruses of public health and veterinary concern. Virus studies in Australian bats usually target the families Paramyxoviridae, Coronaviridae and Rhabdoviridae, with little known about their overall virome composition. We used metatranscriptomic sequencing to characterise the faecal virome of grey-headed flying foxes from three colonies in urban/suburban locations from two Australian states. We identified viruses from three mammalian-infecting (Coronaviridae, Caliciviridae, Retroviridae) and one possible mammalian-infecting (Birnaviridae) family. Of particular interest were a novel bat betacoronavirus (subgenus Nobecovirus) and a novel bat sapovirus (Caliciviridae), the first identified in Australian bats, as well as a potentially exogenous retrovirus. The novel betacoronavirus was detected in two sampling locations 1375 km apart and falls in a viral lineage likely with a long association with bats. This study highlights the utility of unbiased sequencing of faecal samples for identifying novel viruses and revealing broad-scale patterns of virus ecology and evolution.

RNA virome diversity and Wolbachia infection in individual Drosophila simulans flies.

The endosymbiont bacteria of the genus Wolbachia are associated with multiple mutualistic effects on insect biology, including nutritional and antiviral properties. Members of the genus Wolbachia naturally occur in fly species of the genus Drosophila, providing an operational model host for studying how virome composition may be affected by its presence. Drosophila simulans populations can carry a variety of strains of members of the genus Wolbachia, with the wAu strain associated with strong antiviral protection under experimental conditions. We used D. simulans sampled from the Perth Hills, Western Australia, to investigate the potential virus protective effect of the wAu strain of Wolbachia on individual wild-caught flies. Our data revealed no appreciable variation in virus composition and abundance between individuals infected or uninfected with Wolbachia associated with the presence or absence of wAu. However, it remains unclear whether wAu might affect viral infection and host survival by increasing tolerance rather than inducing complete resistance. These data also provide new insights into the natural virome diversity of D. simulans. Despite the small number of individuals sampled, we identified a repertoire of RNA viruses, including nora virus, galbut virus, thika virus and La Jolla virus, that have been identified in other species of the genus Drosophila. Chaq virus-like sequences associated with galbut virus were also detected. In addition, we identified five novel viruses from the families Reoviridae, Tombusviridae, Mitoviridae and Bunyaviridae. Overall, this study highlights the complex interaction between Wolbachia and RNA virus infections and provides a baseline description of the natural virome of D. simulans.

Erratum: Virome composition in marine fish revealed by meta-transcriptomics.

[This corrects the article DOI: 10.1093/ve/veab005.].

Virome composition in marine fish revealed by meta-transcriptomics.

Revealing the determinants of virome composition is central to placing disease emergence in a broader evolutionary context. Fish are the most species-rich group of vertebrates and so provide an ideal model system to study the factors that shape virome compositions and their evolution. We characterized the viromes of nineteen wild-caught species of marine fish using total RNA sequencing (meta-transcriptomics) combined with analyses of sequence and protein structural homology to identify divergent viruses that often evade characterization. From this, we identified twenty-five new vertebrate-associated viruses and a further twenty-two viruses likely associated with fish diet or their microbiomes. The vertebrate-associated viruses identified here included the first fish virus in the Matonaviridae (single-strand, negative-sense RNA virus). Other viruses fell within the Astroviridae, Picornaviridae, Arenaviridae, Reoviridae, Hepadnaviridae, Paramyxoviridae, Rhabdoviridae, Hantaviridae, Filoviridae, and Flaviviridae, and were sometimes phylogenetically distinct from known fish viruses. We also show how key metrics of virome composition-viral richness, abundance, and diversity-can be analysed along with host ecological and biological factors as a means to understand virus ecology. Accordingly, these data suggest that that the vertebrate-associated viromes of the fish sampled here are predominantly shaped by the phylogenetic history (i.e. taxonomic order) of their hosts, along with several biological factors including water temperature, habitat depth, community diversity and swimming behaviour. No such correlations were found for viruses associated with porifera, molluscs, arthropods, fungi, and algae, that are unlikely to replicate in fish hosts. Overall, these data indicate that fish harbour particularly large and complex viromes and the vast majority of fish viromes are undescribed.

Red fox viromes in urban and rural landscapes.

The Red fox (Vulpes vulpes) has established large populations in Australia's urban and rural areas since its introduction following European settlement. The cryptic and highly adaptable nature of foxes allows them to invade cities and live among humans whilst remaining largely unnoticed. Urban living and access to anthropogenic food resources also influence fox ecology. Urban foxes grow larger, live at higher densities, and are more social than their rural counterparts. These ecological changes in urban red foxes are likely to impact the pathogens that they harbour, and foxes could pose a disease risk to humans and other species that share these urban spaces. To investigate this possibility, we used a meta-transcriptomic approach to characterise the virome of urban and rural foxes across the Greater Sydney region in Australia. Urban and rural foxes differed significantly in virome composition, with rural foxes harbouring a greater abundance of viruses compared to their urban counterparts. We identified ten potentially novel vertebrate-associated viruses in both urban and rural foxes, some of which are related to viruses associated with disease in domestic species and humans. These included members of the Astroviridae, Picobirnaviridae, Hepeviridae, and Picornaviridae as well as rabbit haemorrhagic disease virus-2. This study sheds light on the viruses carried by urban and rural foxes and emphasises the need for greater genomic surveillance of foxes and other invasive species at the human-wildlife interface.

Meta-Transcriptomic Identification of Divergent Amnoonviridae in Fish.

Tilapia lake virus (TiLV) has caused mass mortalities in farmed and wild tilapia with serious economic and ecological consequences. Until recently, this virus was the sole member of the Amnoonviridae, a family within the order Articulavirales comprising segmented negative-sense RNA viruses. We sought to identify additional viruses within the Amnoonviridae through total RNA sequencing (meta-transcriptomics) and data mining of published transcriptomes. Accordingly, we sampled marine fish species from both Australia and China and discovered several segments of two new viruses within the Amnoonviridae, tentatively called Flavolineata virus and Piscibus virus, respectively. In addition, by mining vertebrate transcriptome data, we identified nine additional virus transcripts matching to multiple genomic segments of TiLV in both marine and freshwater fish. These new viruses retained sequence conservation with the distantly related Orthomyxoviridae in the RdRp subunit PB1, but formed a distinct and diverse phylogenetic group. These data suggest that the Amnoonviridae have a broad host range within fish and that greater animal sampling will identify additional divergent members of the Articulavirales.

Meta-transcriptomic identification of Trypanosoma spp. in native wildlife species from Australia.

BACKGROUND: Wildlife species carry a remarkable diversity of trypanosomes. The detection of trypanosome infection in native Australian fauna is central to understanding their diversity and host-parasite associations. The implementation of total RNA sequencing (meta-transcriptomics) in trypanosome surveillance and diagnosis provides a powerful methodological approach to better understand the host species distribution of this important group of parasites. METHODS: We implemented a meta-transcriptomic approach to detect trypanosomes in a variety of tissues (brain, liver, lung, skin, gonads) sampled from native Australian wildlife, comprising four marsupials (koala, Phascolarctos cinereus; southern brown bandicoot, Isoodon obesulus; swamp wallaby, Wallabia bicolor; bare-nosed wombat, Vombatus ursinus), one bird (regent honeyeater, Anthochaera phrygia) and one amphibian (eastern dwarf tree frog, Litoria fallax). Samples corresponded to both clinically healthy and diseased individuals. Sequencing reads were de novo assembled into contigs and annotated. The evolutionary relationships among the trypanosomatid sequences identified were determined through phylogenetic analysis of 18S rRNA sequences. RESULTS: We detected trypanosome sequences in all six species of vertebrates sampled, with positive samples in multiple organs and tissues confirmed by PCR. Phylogenetic analysis indicated that the trypanosomes infecting marsupials were related to those previously detected in placental and marsupial mammals, while the trypanosome in the regent honeyeater grouped with avian trypanosomes. In contrast, we provide the first evidence for a trypanosome in the eastern dwarf tree frog that was phylogenetically distinct from those described in other amphibians. CONCLUSIONS: To our knowledge, this is the first meta-transcriptomic analysis of trypanosomes in native Australian wildlife, expanding the known genetic diversity of these important parasites. We demonstrated that RNA sequencing is sufficiently sensitive to detect low numbers of Trypanosoma transcripts and from diverse hosts and tissues types, thereby representing an effective means to detect trypanosomes that are divergent in genome sequence.

A Divergent Articulavirus in an Australian Gecko Identified Using Meta-Transcriptomics and Protein Structure Comparisons.

The discovery of highly divergent RNA viruses is compromised by their limited sequence similarity to known viruses. Evolutionary information obtained from protein structural modelling offers a powerful approach to detect distantly related viruses based on the conservation of tertiary structures in key proteins such as the RNA-dependent RNA polymerase (RdRp). We utilised a template-based approach for protein structure prediction from amino acid sequences to identify distant evolutionary relationships among viruses detected in meta-transcriptomic sequencing data from Australian wildlife. The best predicted protein structural model was compared with the results of similarity searches against protein databases. Using this combination of meta-transcriptomics and protein structure prediction we identified the RdRp (PB1) gene segment of a divergent negative-sense RNA virus, denoted Lauta virus (LTAV), in a native Australian gecko (Gehyra lauta). The presence of this virus was confirmed by PCR and Sanger sequencing. Phylogenetic analysis revealed that Lauta virus likely represents a newly described genus within the family Amnoonviridae, order Articulavirales, that is most closely related to the fish virus Tilapia tilapinevirus (TiLV). These findings provide important insights into the evolution of negative-sense RNA viruses and structural conservation of the viral replicase among members of the order Articulavirales.

Origin, tempo, and mode of the spread of DENV-4 Genotype IIB across the state of São Paulo, Brazil during the 2012-2013 outbreak.

BACKGROUND: Dengue virus type 4 (DENV-4) was first reported in Brazil in 1982 and since then no more cases were detected again in Brazil until 2010, when the virus was reintroduced. Over the following years, the virus spread to several Brazilian states and resulted in about 1,400,000 dengue cases, in 2013. The largest number of cases were documented in the Southeast macro-region. OBJECTIVES: To determine the phylogeography of DENV-4 Genotype IIB strains isolated during the epidemics in 2012-2013 in São Paulo, Brazil, we aimed to contextualise the contribution of viruses sampled in different localities across the overall movement of DENV-4 in Brazil. METHODS: Based on the envelope gene sequences retrieved from GenBank, we employed a Bayesian phylogeographic approach to assess the spatiotemporal dynamics of DENV-4 Genotype IIB in São Paulo, Brazil. FINDINGS: The dispersal dynamics of DENV-4 Genotype IIB in Brazil indicated Rio de Janeiro and Mato Grosso states as the most likely routes toward São Paulo before the 2012-2013 outbreak. Likewise, Guarujá and São José do Rio Preto facilitated viral spread and transmission to other localities in the South and Southeast macro-regions in Brazil. CONCLUSIONS: The spread pattern of DENV-4 Genotype IIB strains across the country supports two independent introductions of the virus in São Paulo in a short period of time. Furthermore, São Paulo appears to have played a pivotal role in the dissemination of DENV-4 to other locations in Brazil.

Origin of the São Paulo Yellow Fever epidemic of 2017-2018 revealed through molecular epidemiological analysis of fatal cases.

The largest outbreak of yellow fever of the 21st century in the Americas began in 2016, with intense circulation in the southeastern states of Brazil, particularly in sylvatic environments near densely populated areas including the metropolitan region of São Paulo city (MRSP) during 2017-2018. Herein, we describe the origin and molecular epidemiology of yellow fever virus (YFV) during this outbreak inferred from 36 full genome sequences taken from individuals who died following infection with zoonotic YFV. Our analysis revealed that these deaths were due to three genetic variants of sylvatic YFV that belong the South American I genotype and that were related to viruses previously isolated in 2017 from other locations in Brazil (Minas Gerais, Espírito Santo, Bahia and Rio de Janeiro states). Each variant represented an independent virus introduction into the MRSP. Phylogeographic and geopositioning analyses suggested that the virus moved around the peri-urban area without detectable human-to-human transmission, and towards the Atlantic rain forest causing human spill-over in nearby cities, yet in the absence of sustained viral transmission in the urban environment.

Origin, tempo, and mode of the spread of DENV-4 Genotype IIB across the state of São Paulo, Brazil during the 2012-2013 outbreak

BACKGROUND Dengue virus type 4 (DENV-4) was first reported in Brazil in 1982 and since then no more cases were detected again in Brazil until 2010, when the virus was reintroduced. Over the following years, the virus spread to several Brazilian states and resulted in about 1,400,000 dengue cases, in 2013. The largest number of cases were documented in the Southeast macro-region. OBJECTIVES To determine the phylogeography of DENV-4 Genotype IIB strains isolated during the epidemics in 2012-2013 in São Paulo, Brazil, we aimed to contextualise the contribution of viruses sampled in different localities across the overall movement of DENV-4 in Brazil. METHODS Based on the envelope gene sequences retrieved from GenBank, we employed a Bayesian phylogeographic approach to assess the spatiotemporal dynamics of DENV-4 Genotype IIB in São Paulo, Brazil. FINDINGS The dispersal dynamics of DENV-4 Genotype IIB in Brazil indicated Rio de Janeiro and Mato Grosso states as the most likely routes toward São Paulo before the 2012-2013 outbreak. Likewise, Guarujá and São José do Rio Preto facilitated viral spread and transmission to other localities in the South and Southeast macro-regions in Brazil. CONCLUSIONS The spread pattern of DENV-4 Genotype IIB strains across the country supports two independent introductions of the virus in São Paulo in a short period of time. Furthermore, São Paulo appears to have played a pivotal role in the dissemination of DENV-4 to other locations in Brazil.

Evolutionary history and spatio-temporal dynamics of dengue virus serotypes in an endemic region of Colombia.

Dengue is a prevalent disease in Colombia and all dengue virus serotypes (DENV-1 to -4) co-circulate in the country since 2001. However, the relative impact of gene flow and local diversification on epidemic dynamics is unknown due to heterogeneous sampling and lack of sufficient genetic data. The region of Santander is one of the areas with the highest incidence of dengue in Colombia. To provide a better understanding of the epidemiology of dengue, we inferred DENV population dynamics using samples collected between 1998 and 2015. We used Bayesian phylogenetic analysis and included 143 new envelope gene sequences from Colombia, mainly from the region of Santander, and 235 published sequences from representative countries in the Americas. We documented one single genotype for each serotype but multiple introductions. Whereas the majority of DENV-1, DENV-2, and DENV-4 strains fell into one single lineage, DENV-3 strains fell into two distinct lineages that co-circulated. The inferred times to the most recent common ancestors for the most recent clades of DENV-1, DENV-2, and DENV-4 fell between 1977 and 1987, and for DENV-3 was around 1995. Demographic reconstructions suggested a gradual increase in viral diversity over time. A phylogeographical analysis underscored that Colombia mainly receives viral lineages and a significant diffusion route between Colombia and Venezuela. Our findings contribute to a better understanding of the viral diversity and dengue epidemiology in Colombia.

Codon adaptation biases among sylvatic and urban genotypes of Dengue virus type 2.

Dengue virus (DENV) emerged from the sylvatic environment and colonized urban settings, being sustained in a human-Aedes-human transmission chain, mainly by the bites of females of the anthropophilic species Aedes aegypti. Herein, we sought evidence for fine-tuning in viral codon usage, possibly due to viral adaptation to human transmission. We compared the codon adaptation of DENV serotype 2 (DENV-2) genotypes from urban and sylvatic habitats and tried to correlate the findings with key evolutionary determinants. We found that DENV-2 codons of urban and sylvatic genotypes had a higher CAI to humans than to Ae. aegypti. Remarkably, we found no significant differences in codon adaptation to human between urban American/Asian and sylvatic DENV-2 genotypes. Moreover, CAI values were significantly different, when comparing all genotypes to Ae. aegypti codon preferences, with lower values for sylvatic than urban genotypes. In summary, our findings suggest the presence of a molecular signature among the genotypes that circulate in sylvatic and urban environments, and may help explain the trafficking of DENV-2 strains to an urban cycle.

First Genome Sequences of Dengue Virus (DENV) Strains Isolated during the First DENV-4 Outbreak in São Paulo, Brazil.

Dengue virus (DENV) is an arbovirus belonging to the genus Flavivirus, family Flaviviridae In Brazil, the reemergence and spread of DENV type 4 (DENV-4) across the country were responsible for a significant outbreak in Guarujá, São Paulo, Brazil. Here, we report the first genomic sequences of DENV strains circulating in Guarujá during the 2013 outbreak.

Virus del dengue de serotipo 1 (DENV-1) de Colombia: su contribución a la presentación del dengue en el departamento de Santander / Dengue virus serotype 1 (DENV-1) from Colombia: its contribution to dengue occurrence in Santander

Introducción. Los cuatro serotipos del virus del dengue circularon en el departamento de Santander entre 1998 y 2008. No existe información sobre el papel del serotipo 1 (DENV-1) en la epidemiología de la enfermedad. Objetivo. Analizar la relación entre el cambio de predominancia del (DENV-1) con su diversificación genética, predominancia de los otros serotipos y presentación del dengue grave. Materiales y métodos. La diversificación genética se estudió por análisis filogenético usando la secuencia del gen E de 12 cepas del virus. Para el análisis se utilizaron datos sobre predominancia de los serotipos obtenidos en estudios previos y datos oficiales de incidencia del dengue. Resultados. Los virus seleccionados se agruparon en el genotipo V junto a (DENV-1) de países de Latinoamérica y se evidenció segregación en cuatro linajes. Los cambios en la predominancia del virus coincidieron con el reemplazo de linaje y esto, a su vez, con incremento en la prevalencia de DENV-2 y DENV-3, e incremento del dengue grave. Conclusión. La diversificación genética podría contribuir a cambios de predominancia de (DENV-1), y la relación del virus con el DENV-2 y DENV-3 en situaciones que favorecen la presentación de casos graves. Se necesitan más estudios para precisar el papel de los serotipos en la epidemiología del dengue.

Introduction: Between 1998 and 2008 all dengue virus serotypes circulated in the Departamento de Santander, an endemic region in northeastern Colombia. No information is available as to the role of serotype 1 (DENV-1) with respect to epidemiology of dengue. Objective: To analyze the relationship between changes in DENV-1 predominance with respect to genetic diversity, prevalence of others serotypes and occurrence of severe dengue. Methods: Virus genetic diversity was studied by phylogenetic analysis comparing E gene sequences from 12 viral strains. Data about serotypes predominance obtained in previous studies and official data about dengue incidence were used for analysis. Results: Selected viruses grouped into genotype V together DENV-1 from Latin America countries, and segregation in four lineages was evidenced. Changes in virus predominance coincided with replacement of lineage, increase in prevalence of DENV-2 and DENV-3 and increase of severe dengue. Conclusion: Genetic divergence could have contributed to changes in DENV-1 predominance. The relationship of the virus with DENV-2 and DENV-3 could create scenarios that promote occurrence of severe cases. More studies are required to ascertain the precise role of serotypes in the epidemiology of dengue.

[Dengue virus serotype 1 (DENV-1) from Colombia: its contribution to dengue occurrence in Santander]. / Virus del dengue de serotipo 1 (DENV-1) de Colombia: su contribución a la presentación del dengue en el departamento de Santander.

INTRODUCTION: Between 1998 and 2008 all dengue virus serotypes circulated in the Departamento de Santander, an endemic region in northeastern Colombia. No information is available as to the role of serotype 1 (DENV-1) with respect to epidemiology of dengue. OBJECTIVE: To analyze the relationship between changes in DENV-1 predominance with respect to genetic diversity, prevalence of others serotypes and occurrence of severe dengue. METHODS: Virus genetic diversity was studied by phylogenetic analysis comparing E gene sequences from 12 viral strains. Data about serotypes predominance obtained in previous studies and official data about dengue incidence were used for analysis. RESULTS: Selected viruses grouped into genotype V together DENV-1 from Latin America countries, and segregation in four lineages was evidenced. Changes in virus predominance coincided with replacement of lineage, increase in prevalence of DENV-2 and DENV-3 and increase of severe dengue. CONCLUSION: Genetic divergence could have contributed to changes in DENV-1 predominance. The relationship of the virus with DENV-2 and DENV-3 could create scenarios that promote occurrence of severe cases. More studies are required to ascertain the precise role of serotypes in the epidemiology of dengue.